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Abstract 



We construct lepton flavour-dependent and basis-independent measures of i?-parity 
violation in the minimal supersymmetric extension of the Standard Model (MSSM). We 
Qh. work in the context of exact super symmetry, neglecting the effects of Higgs vacuum ex- 

pectation values and soft super symmetry-breaking terms. We devote particular attention 
to appropriate choices of flavour eigenstates, and to the counting and enumeration of 
i?-violating invariants in two- and three-generation models. We also make an indicative 
application of our results to derive possible basis-independent cosmological upper bounds 



on flavour-dependent violation of R parity. 
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1 Introduction 

One of the most attractive possible extensions of the Standard Model is super symmetry (SUSY) 
[|TJ. The Lagrangian for the minimal supersymmetric extension of the Standard Model (MSSM) 
could in principle contain renormalisable baryon- (B-) and lepton- (L-) number-violating in- 
teractions. These are constrained by the fact that neither B nor L violation has yet been ob- 
served experimentally. It is possible to remove these interactions by imposing a multiplicatively- 
conserved global symmetry called R parity |2|], whose action can be defined at the level of the 
MSSM superpotential such that matter superfields (L,Q,E C ,U C ,D C ) are multiplied by — 1, 
but the Higgs superfields are invariant. In component notation, R parity multiplies fields by 
(— 1) 3B+L+2SI , where S is the spin. 

Alternatively, one can merely require that possible R- violating couplings be sufficiently small 
that their effects would not (yet) have been seen by experiment 0. One might also want to 
bear in mind the drastic cosmological implications of possible R violation, even at levels below 
the current experimental upper limits, for example for the survival of any primordial baryon 
or lepton asymmetry |3], ||. However, there is problem with the formulation of experimental 
or cosmological upper bounds on the magnitudes of R- violating couplings: L violation can be 
moved around the superpotential by basis transformations, so the apparent magnitudes of these 
coupling constants are basis dependent. In a previous paper || we have therefore derived basis- 
independent expressions that parametrize R violation in the MSSM: we denote || subsequently 
asDE. 

The /^-violating invariants of DE have various limitations: for example, there we neglected 
soft SUSY-breaking interactions and the Higgs vacuum expectation values, which we do here 
also. Another shortcoming of DE was that we summed over lepton flavours, since otherwise the 
invariants would have depended on the basis chosen for the singlet right-handed leptons. This 
means that the basis-independent expressions of DE only measured the sum of the violation 
of the lepton flavours L e , L^ and L T . In this paper, we determine the appropriate basis for 
the quarks and singlet right-handed leptons, obviating the need to sum over the quark and 
singlet lepton flavours. In this way, we obtain basis-independent measures of the violation of 
individual lepton flavours. Since there are more terms in the flavour sums of DE than there are 
R- violating couplings, not all the terms in our previous sums can be independent invariants. We 
therefore discuss the number of distinct, basis-independent combinations of coupling constants 
that measure i?-parity breaking if one does not sum over lepton flavours, and express the 
"superfluous" invariants in terms of an independent basis set. We then use these to extend our 
previous discussion of the possible cosmological constraints on i?-parity violation to include the 
case where additional flavour symmetries are present. 

2 Review of Flavour- and Basis-Independent Measures 
of R Violation 

The problem that DE sought to address is that the magnitude of possible R- violating superpo- 
tential couplings is basis dependent. The MSSM superpotential, including possible .R-violating 



terms, is: 



W = ^ H HH+h^HQ i U^+h^HQ i D c 1 +h^HL i E^+e i HL i +\fL i L j E c k +\fL i Q 3 Dl+\fU^D c j D c k 

(1) 
Hereon, we focus on the violation of lepton number, and so neglect the possible 5-violating 

coupling constant A3. It is well known that e in ([!]) can be rotated away by a basis transfor- 
mation. This is because the doublet Higgs and lepton superfields H, L t have the same gauge 
quantum numbers, and SUSY removes the Standard Model distinction between H and the L, 
that is based on their spins. Once we allow for R violation, there is no unambiguous distinction 
between H and the Li. If it is unclear how to define a lepton, then it is also unclear which 
coupling constants violate L. 

If the Higgs and doublet lepton superfields are assembled into a vector 

</>i={H,Li) (2) 

whose index I runs from to 3, then the superpotential ([[]) can be rewritten as: 

W = n'Hfa + Xi Jr MjE c k + X'rhQpDl + KHQJJ'r (3) 

where /i 7 = (//,£*), \h{ k = X°J k , Xf = X^ k , X° d jk = hf, and Xf = Xf. This superpotential 
contains R violation if different coupling constants, as vectors and matrices in space, choose 
different directions to be the Higgs. The differences between the various directions chosen to 
be the Higgs parametrize the amount of R violation present, and are independent of basis. 
Therefore, to obtain basis- independent measures of R violation ||, one should first consider a 
lepton Yukawa interaction, which defines a plane spanned by a Higgs and lepton in <j) space, 
then use another interaction to fix one direction in that plane as the Higgs, with the orthogonal 
direction being defined as the lepton. If some other interaction chooses as the Higgs a direction 
that has a non-zero projection onto the direction previously chosen to be a lepton, this yields 
a measure of R violation that is independent of the basis chosen in space. Note that at least 
one Yukawa coupling must enter into any such measure of R violation, because, in the absence 
of the Yukawa couplings, all the directions in space are equally good choices to be the Higges, 
and it is not possible to define R violation. 

A number of algebraic expressions, which correspond to sums of closed supergraphs, were 
computed in DE using this algorithm: 

{^xj J n J d n^ K ^ Lp X qr ) 

|/i| 2 |A e | 2 |A 

/i J *A^PA^Af Lg A^ M P> M -l/2(/i 7 *A^PA e J ^> K )(A^ M ^A^) 

\^Tr[XTXlXTXV l)) 



6l = I..I2IX 1211.12 ( 4 ) 



\e 



\Irs* \IJp \JKq* \KLq \LMp* \Mrs i In/ \Irs* \IJp\JKq* \Krs\( \LMp* \LMq\ 

\X d \>Tr[XTX P eXfxl] " ' (6) 

Tr[XlXTXlXl*XlXT] l ' 

1 Tr[AgAf X e \?]Tr[\%K*} + Tr[A^Af AgAr]Tr-[AgAf] 

2 Tr[AfAfA-ArAlAf] 



* 5 = wW • (8) 

Here the traces are over the capitalized (f) indices, and the lower-case indices correspond to 
quark and right-handed lepton generations, which are also summed []. 

Some of the couplings in eqn (|3|) yield interactions of particles in the subspace with 
particles in other multi-dimensional spaces, e.g. E c , Q, D c . The magnitudes of the couplings 
therefore depend on the choice of basis in these other spaces. In DE, this was resolved by 
summing over all the indices: those of 0, of the quarks, and of the right-handed leptons. 
This produced basis-independent expressions that measured the total L violation. However, 
expressions that measure the magnitude of R violation for a given quark or lepton generation 
would have additional phenomenological applications. It would also be interesting to know 
how many independent invariants there are: it is clear that, if one did not do the quark and 
right-handed lepton sums in the expressions (|J) -@, there would be more "invariants" than 
R- violating couplings, so some of these invariants must be redundant. 

3 Specification of the Superfield Bases 

There are various possible choices for the Q, D c and E c bases. The most useful ones are 
presumably those that diagonalise the propagators in flavour space. In the absence of soft 
SUSY breaking, this diagonalisation can be performed at the superfield level. Moreover, we are 
neglecting the Higgs vacuum expectation values, as is appropriate for conditions in the early 
Universe. We may therefore choose the Q basis such that 

E * W p * + £ hi'hzr « s« r , (9) 

Ip a 

EAfAfocr, (10) 

i P 

£ \l Jg \i Jr * oc 5 qr . (11) 

IJ 

These conditions ensure that the sums of the vacuum polarization supergraphs are flavour 
diagonal for each of the corresponding types of particle:Q, D c , E c . For example, the conditions 
(|S], [IT], [TTJ) ensure that the one-loop propagators are flavour-diagonal at high temperatures in 
the early universe, when soft SUSY breaking and Higgs vacuum expectation values can be 
neglected. In the absence of -R-parity breaking, the above prescription yields the Standard 
Model mass-eigenstate basis for the E c and D c : writing X°J k = \h? k (the Yukawa coupling 
matrix of the leptons) and X l J k = 0, equation flllf) is equivalent to h\h e being diagonal. This 

^^Note that we have rewritten slightly 84, as compared to the expression in DE. If one sums over all the 
right-handed lepton indices, the expression above is equivalent that in to DE. However, if one does not perform 
the right-handed index sum, the above expression vanishes in the absence of R violation, whereas the version 
in DE does not. 



the D c basis such that 



and the E c basis such that 



is the usual mass-eigenstate basis for the singlet right-handed charged leptons in the Standard 
Model, and a similar argument applies to D c . 

4 Counting the Flavour-Dependent and Basis-Independent 
Invariants in Two Generations 

Having chosen a well-motivated basis for the quarks and right-handed leptons, we can now 
construct invariants that measure the violation of individual lepton flavours, where these are 
defined to be associated with the right-handed leptons E^. To illustrate how this may be done, 
we first consider as a warm-up exercise a simplified model with just two generations, whose 
geometry is easier to visualise than the fully realistic three-generation case. If one simply 
considers all the expressions listed in equations (f|) - (|j), and does not perform the quark and 
singlet-lepton sums, one obtains 42 invariants in the two-generation model: eight from Si, two 
from 82, eight from £3, none from £4, and 24 from 65 0. It is clear that this is overcounting, 
because in two generations there are only twelve possible .R-violating couplings: two from each 
of \x and the four A^, and one from each A e Q. The reason that there are more invariants 
than R- violating couplings is as follows. The different couplings /i, A^, A e , ... choose different 
directions in <fi space to be the Higgs and the leptons, as seen in Fig. [I]. The various invariants 
measure the differences between these vectors, which is a geometric invariant independent of 
the basis chosen in the <p space. However, it is geometrically clear that if one knows the vectors 
in space corresponding to the differences between \i and A^ 1 and between \x and A^ 2 , then the 
invariant measuring the difference between A^ 1 and A^ 2 is not providing any new information. 

Given the excess of invariants, we now analyze the number of independent invariants, and 
express the extra dependent ones in terms of an independent basis set. This is simpler in two 
generations, because the space is three-dimensional, and hence easier to visualize. We expect 
to find 10 independent invariants in two generations, because, if one chooses to work in the 
basis where fi = (fi, 0,0), two L- violating couplings arise from each of the four X^ q , and one 
L- violating coupling from each of the A e . This expectation of finding ten independent invariants 
is indeed confirmed by our analysis. 

The A e have antisymmetric indices, so they have a zero eigenvector in this two-dimensional 
model. This means that they can be pictured as a plane in space, spanned by the Higgs and 
a lepton: 

\l Jk = IX'KH 1 L( - LlH J ) (12) 

where |A*| = Tr[A*A^]. The intersection of the two planes chooses a direction for the Higgs, 
which is here labelled H. The directions orthogonal to H in the planes defined by the \ l e are 
labelled by Li (i = 1, 2). The lepton directions Li are orthogonal by eqn (O): 



= A^ J1 Af 2 oc (H'Ll - l[E J ){E J l{ - L J 2 H r ) = -2|#|% • L 2 (13) 

These definitions of H and Li provide the axes in Fig. [I]. 

2 The corresponding number in the full three-generation model would be 321 invariants. 
3 There are 39 such couplings in the three-generation model. 




Figure 1: Directions chosen by the interactions //, X^ q and A* in the two-generation space: 
i,p,q = 1,2. The lepton Yukawa couplings A* and \ 2 e each choose a plane which is spanned 
by the Higgs and a lepton. The A e choice of the Higgs is therefore the intersection of the two 
planes, here labelled H. In general, \x and the X^ q each choose a preferred direction to identify as 
the Higgs: their projections onto the choice of L\ and L 2 made by the A e are basis-independent 
measures of .R-parity violation. Note that there are in principle four distinct X^ Q : for simplicity, 
we show only one. 



There is R violation in this two-generation model if the direction chosen to be the Higgs by 
the A e couplings differs from the direction chosen by jjl or the A^ couplings. The lepton-number 
violation for each flavour can be parametrized by taking the inner product between \x or A^ - 
directions in <fi space that should correspond to the Higgs - and Li as defined by the A e . This 
gives the following 10 independent invariants 
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(17) 
where the first expression in each case indicates the geometric interpretation of the invariant as 
indicated in Fig. [1], and the A e A e traces in the second terms in all of these formulae are zero in 
the basis chosen by (|TT|) . The subscripts are derived from the numbering of the 5 invariants in 
equations (f|) to ([8]), and the superscripts are the generation indices of the right-handed leptons 
and quarks. One can also construct other invariants, corresponding to 8\ and 85 without their 
generation sums. However, these can be expressed in terms of the invariants (|T^) - flT7|), as 
shown below, and so are not independent. We therefore find the expected 10 invariant measures 
of R violation in this two-generation model. We can express other possible invariants in terms 
of the independent basis set ( |T^D -- (|T7|). The combination of coupling constants that appears 
in the relation between one set of invariants and the other is y8, which may be interpreted as 
the length of a vector in space, see (|14]) - - fjl7[) . For instance, the vector corresponding to 
the difference between X p d and X r d s , is the vector sum of the differences between X d q and /i, and 
between \i and X r d s . 

In the basis chosen by the A e , the invariants can written as: 

(A£Vo-^(Af)o 
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It is clear from fll8|) that one can express 



8? q = V^ M (1 - % ~ #) ~ V^(l - %** - % m ) ,j¥>i (22) 



and from (]2lf) one finds 



cpqirs 



$p«(i - ST a - 5i irs ) 
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J ^i 



(23) 



where we have written /io = |/i|[l — 8$ — Sl 1 } 1 ^ 2 , for instance. 

5 Extension to the Three-Generation Case 

The full three-generation model is more complicated to analyze, both because it is more difficult 
to visualize - the <fi space is four-dimensional, and because the A e can contain R violation 
independently of the other interactions. We therefore count invariants by making a particular 
basis choice that is motivated by the coupling constants. The couplings /i and A^ are still vectors 
in space that would like to define the direction corresponding to the Higgs, but each A e is now 
geometrically represented by two non-intersecting planes. To see this, we recall that, in three 
generations, any A e (say A*) is a four-dimensional anti-symmetric matrix whose eigenvalues can 
be labelled as ±ihi and ±ic\. This means that there is a basis in which 
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(24) 



Each of the two non-intersecting planes in space represented by ( |24D would like to be spanned 
by a Higgs and a lepton. However, if a pair of the Ag eigenvalues vanish, then A e corresponds 
to a single plane, as in the two-generation model discussed in the previous section, and lepton 
flavour is conserved. In the general non-degenerate case, both of these planes should be spanned 
by a Higgs and a lepton field. This means that even a single A e coupling does not in general 
conserve lepton flavour, because one right-handed lepton can interact with two flavours of 
left-handed lepton. Moreover, it would also like two directions in space to be the Higgs. 

In a similar fashion, a second coupling Ag can also be pictured as two planes. The directions 
in which the Ag planes intersect with the Ag planes can be chosen as the Higgs and, for instance, 
L 3 . The other leptons L x and L 2 can be taken as the orthogonal directions in the Ag planes. 
With this choice, the Higgs and L 3 fields do not interact, and nor do L\ and L 2 , so Ag takes 
the form 

ei e 2 

-ei -6a 

-e 2 6i 
b 2 -&i 



A 



(25) 



The third coupling of this type, Ag, is of generic form in the basis chosen in this way by Ag and 
Ag. The counting of lepton-flavour- violating couplings therefore yields one from Ag, three from 
Ag and five from Ag. However, there are also three conditions imposed on these couplings by 
the choice of right-handed lepton basis (|TT|) , so we expect 6 independent invariant measures of 
R violation 



We could take, for instance, 51 



Ajk 



1, 2, 3 and <5 4 ,i ^ j ^ k. This choice has 



the feature that the Sf l parametrize Li violation in the A*, and the remaining three invariants 
then measure the differences between the bases chosen by the various A*. 

As in the two-generation case, /i and the X p d q choose directions in <f> space that should be 
identified with the Higgs. The components of these vectors along the directions chosen by 
the A e to be leptons are invariant measures of R violation. We therefore need three further 
invariants for each of /i and the X d q , which can be taken to be: 
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(26) 

(27) 
(28) 

(29) 
(30) 
(31) 



Note that |/zAgA 3 | 2 is not included in the above list, although |/zAgA 2 | 2 is. This is because the 
role of the expressions (^) - ([!!]) is to determine whether fi and/or A^ share the same definition 
of the Higgs with the A e , whereas the differences among the A e have already been accounted 
for by the <5 4 invariants. The purpose of the second A e factor in expressions ( [Z6"D - ([JI]) is only 
to choose a direction in the plane of the first A e to be the Higgs. The orthogonal direction is 
then the lepton, and one projects \i or A^ onto the direction of this lepton field. If the three A e 
happen to conserve R parity, which occurs if they each represent a single plane, and the three 
planes intersect in a common direction, which is then defined unambiguously to be the Higgs, 
then /iAgA 3 ~ fiihih 3 and yuA^A 2 ~ fiihih 2 . Thus one measures the same i?-violating coupling, 
whether the second matrix is taken to be X 2 or A 3 . We therefore have 6 + 30 = 36 invariants 
in the three-generation case, as expected. 

Another way to count invariants in the three-generation case, which may seem more trans- 
parent, is to choose the Higgs to be the direction corresponding to /x, and the leptons to be 
Li = n ■ A*. The right-handed lepton basis must then be chosen such that 



/iU^A^/i oc 5 lJ 
so that the Lj are orthogonal in space. The A e are then of the form 
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The 9 x 3 = 27 invariants 
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then measure Li violation in the X^ q , and the nine invariants 
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measure Lj violation in X l e : 
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(36) 



completing our second enumeration of a basis of independent flavour-dependent measures of R 
violation 0. 

6 Application to Lepton Flavour Violation in Early Cos- 
mology 

In the hot plasma of the early Universe at temperatures above the electroweak phase transition 
(EPT): T zL T c ~ 100 GeV, non-perturbative processes that violate 5 + i are in thermal 
equilibrium. If i?-parity non-conserving interactions that violate all the Lj are also in thermal 
equilibrium above the EPT, then any asymmetry in B or the Lj gets washed out. This means 
that either the baryon asymmetry of the Universe (BAU) observed today was not present in the 
thermal plasma above the electroweak phase transition, motivating scenarios in which it was 
made at the EPT, or the i?-violating couplings are sufficiently small to be out of equilibrium 
above the EPT. This argument has been used to place strong constraints on .R-violating Yukawa 
couplings ||, |5j of order 

A K <10- 7 (37) 

where \r is some generic i?-violating Yukawa coupling. In this section we review these con- 
straints in the basis-independent approach of DE, with the added ingredient that we discuss 



4 Note that, in the two-generation case, A| can not violate L and Lj simultaneously. This is because \\ 
couples E% to either HL\, HL<z or L\Li- Therefore, in the two-generation case it is easy to check that S^ 3 = 0. 
However, in the three-generation case, X l J k can violate L and L^ simultaneously via the couplings X l J k i,j, ^ k, 
and this is measured by S^ 3 . 
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carefully the bounds on each lepton flavour separately, using the flavour- dependent invariants 
isolated in the previous sections (]. 

We first expand on the previous paragraph's discussion of the origin of the bounds. The 
minimal supersymmetric extension of the Standard Model (MSSM) has four global quantum 
numbers: B,Li,L2,Lz that are conserved at the perturbative level, in the absence of the R- 
violating couplings ([!]). However, anomalous non-perturbative processes that are in thermal 
equilibrium above the electroweak phase transition, but irrelevantly weak today, violate the 
combination B + L, so the global conserved quantum numbers in the MSSM are B/3 — Lf 
i : 1,2,3. If one then introduces the R- violating interactions of ([!]), there are no conserved 
global quantum numbers left. This means that, with the MSSM particle content, all asymme- 
tries in the plasma are washed out before the electroweak phase transition if the i?-violating 
interactions are in thermal equilibrium. As has already been mentioned, one possible origin 
for the baryon asymmetry observed in the Universe today, is to make the asymmetry at the 
transition, which may be possible for certain areas of SUSY parameter space HI. One can also 



create the asymmetry from particles decaying after the EPT || [K| , or store it in particles not 

in thermal equilibrium with the MSSM while the i?-parity violating interactions are present 

12^1 . For the rest of this paper, we will neglect these scenarios, and assume that the BAU was 



present in the plasma above the EPT, and must be protected from the depredations of the 
i?-violating interactions. 

We assume that non-perturbative B + L violation is in thermal equilibrium before the EPT, 
and consider only the renormalisable .R-violating interactions of equation (HD 0. A sufficient 
condition for the BAU to survive is that all the .R-parity violating interactions should be out 
of thermal equilibrium before the EPT, which gives the bound fl37l) . However, this is too 
restrictive, because it means that the asymmetries in all three of the B/3 — Li are preserved, 
whereas keeping one of them would be sufficient to preserve the BAU. 

More precise bounds can be obtained by using our invariants in an analysis based on chemical 
potentials Q We assume an interaction to be in effective chemical equilibrium in the early 
Universe if the interaction rate exceeds the expansion rate. If this condition is satisfied, the 
sum of the chemical potentials of the particles participating in the interaction is zero [jnj| 0. In 



the Standard Model and the MSSM, all the interactions are in equilibrium once the temperature 



has dropped to ~ 1 — 10 TeV, imposing |14| a set of linear equations on the chemical potentials 
that we now review. Since the gauge bosons have zero chemical potential, all the members of 
a gauge multiplet share the same chemical potential. The Higgs interactions in the Standard 
Model imply: 

(38) 
(39) 
(40) 



q l 


+ h- 


-u 3 R 


= 


Q 1 


-h- 


-d? R 


= 


e 


-h- 


-e* 


= 



5 For a more extensive discussion in the mass-eigenstate basis, see ||. 

6 Non-renormalisable interactions are different [fL3|], in that they may be out of equilibrium before all the 
Standard Model interactions come into equilibrium, in which case the asymmetries are not neccessarily washed 
out. 

We follow her the naive kinetic theory discussion of |q, |l4f| , which is adequate for our purpose. For a more 
complete analysis using the grand canonical ensemble, see pd] . 

8 We work in a mass-eigenstate basis. 
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where the chemical potential of a particle is labelled by the particle's name: q is the left-handed 
quark doublet, h the Higgs doublet, and so on. The flavour non-diagonal quark couplings to the 
Higgs imply that all the left-handed quarks have the same chemical potential, so q l = g J = q, 
and similarly for the different flavours of ur and dft quarks. The anomalous B + L violation 
transforms left-handed quarks into left-handed anti-leptons: 

3iV 5 g + ]Tf = (41) 

i 

where N g is the number of quark generations. We therefore have six linear equations for 10 
unknowns. The four free parameters correspond to the four conserved charges in the plasma: 
Qem and the B/3—Li. The plasma can contain asymmetries associated with the three conserved 
global quantum numbers B/3 — Li, but is required not to have an electric charge f\ asymmetry. 
The electric charge density can be written as 



ii, 



Y,Qa{n a -n- a ) (42) 



where the sum over a runs over all the particles in the plasma, Q a is the particle electric 
charge, and n a (n a ) is the (anti-)particle number density. For masses and chemical potentials 
much smaller than the temperature, the particle asymmetry is : 



- 9^ 

n — n ~ 



2 s 

m 



1 + \p 



and 

9^ 



n — n 



m 



i + o h* 



for bosons, (43) 



for fermions, (44) 



6 

where g is the number of internal degrees of freedom of the particle, and m is the zero- 
temperature mass 0. If we neglect the mass contributions in (|3) and (0), since their contri- 



bution to the asymmetry is reduced with respect to the first terms by m 2 /T 2 <C 1, the electric 
charge density in the Standard Model is 



n<- 



N g (q + 2u R -d R )- y £,£ i -J2e) t + 2h = (45) 



(This equation must be multiplied by three in the MSSM at T > m S usY-) In the MSSM, 
there are more than twice as many new particles, and therefore many new chemical potentials, 
but these turn out all to be equal to those of the Standard Model particles. The gauginos are 
majorana particles, so can not store any asymmetry and must have zero chemical potential. 
The gaugino-scalar-fermion interactions then imply that the chemical potentials of the spartners 
(= 4>) are equal to those of the Standard Model particles(= ip): 

<p = ip (46) 



9 Or hypercharge - it amounts to the same constraint |l4j. 

10 See |1(J ||, [16| for an analysis including these mass effects, and [[llj for an careful discussion of the relation 
of the chemical potentials and masses to the asymmetries. We use here [i as a generic chemical potential, and 
hope that the distinction between the chemical potential and the /x term in the MSSM superpotential will be 
clear from the context. 
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The only remaining new chemical potential is that of the second Higgs H, which is required by 
the mass term fiHH to have equal and opposite chemical potential to the Higgs H. 

The solution of the equations of chemical equilibrium in the absence of R violation is eas- 
ily obtained. One can represent the chemical potentials in the MSSM as a vector in a ten- 
dimensional space: ft = (i 1 ,^, £ 3 , e R , e 2 R , e 3 R , h, q, ur, da). The equilibrium conditions (|38| - f|l 



and (|45|) force this vector to be orthogonal to the directions associated with Q em and the in- 
teractions: for instance, the interaction (|38|) corresponds to the direction (0,0,0,0,0,0,1,1,-1,0). 
The conservation of the three B/3 — L t implies that the projections of the chemical potential 
vector fl onto the directions corresponding to B/3 — Li are constant, fixing fl. The direction 
corresponding to B/3 — L\ is 

i7|_ Ll = (-2,0,0,-l,0,0,2,l,l) (47) 

and there are similar expressions for vb_ L2 and vb_ l . it is important to note that these vectors 



are not orthogonal to the interaction vectors fl38[ - |41|) , nor to the electric charge direction fl45|). 
This is to be expected, since the Yukawa interactions can change the ratio of left- and right- 
handed quarks, although they cannot change the total quark number. One cannot simply 
express fl in terms of the vb_ l _. Instead, one can subtract from the vb_ l . the components 



along the directions corresponding to (^ - [11]) and ([15|), leaving three vectors V{ that are the 



projections of the vb_ l , onto the space allowed for fl. Then fl can be expressed in terms of the 

Vi, as: 

n=riB Li Vi (48) 

3 l 

where tjb_ l , is the conserved ratio of B/3 — L, to photons. The baryon-to-photon ratio is 



therefore 1 1 



Vb = fl- B = Y,Vf- L V i - B = ^Y,Vf- Ll (49) 

where B = (0, 0, 0, 0, 0, 0, 0, 6, 3, 3). 

Now suppose we include the i?-violating interaction A3 UfD'jD'f, of equation ([!]) to the 
thermal soup above the electroweak phase transition. We assume for the moment that all the 
Lj-violating couplings are zero or small enough to satisfy (p7|), and neglect them. If A3 is in 
thermal equilibrium, i.e., at least one Xijk > 10~ 7 in the thermal mass-eigenstate basis, then 
there is an additional constraint on the chemical potentials 

ur + d R + d R = . (50) 

In conjunction with (^ - f|l|) and (^), this constraint forces q = 0, so there is no baryon 
asymmetry. An alternative way to formulate this is to note that the direction in chemical 
potential space corresponding to B can be written as a linear combination of the directions (^ 
- 0), fl50"D and ([45]), to which fl must be orthogonal. So, in the presence of the non-perturbative 



B + L violation and even one of the A^ , the baryon asymmetry will be washed out, up to lepton 
mass effects [|T0 , |i~6[| 0. The baryon-number-violating couplings A3 3 should therefore satisfy 



n This caveat applies when there are different asymmetries in the various lepton flavours, but is unlikely to 

13 



This is unsurprising, since any one of the A3 interactions would violate all the B/3 — Li, 
at least one of which needs to be conserved Q. 

We now consider the Lj-violating couplings present in ([[]), assuming that the A3 3 satisfy 
(pTFp. We wish to find the minimal bounds on the various lepton-flavour- violating rates that 
ensure that an asymmetry in at least one of the B/3 — Li is preserved. If we include interactions 
that violate all of L\, L2 and L3, they imply the chemical potential relations 

h = e (51) 

for i = 1,2, 3, implying that all the asymmetries are washed out. Even the lepton mass effects 
cannot help in this case, because all the lepton flavours would have the same asymmetry by 
fl5"l"|) . However, if the interactions in equilibrium conserve a lepton flavour, say L\ = L e , then 
the asymmetry in B/3 — L\ is conserved, and a baryon asymmetry 

28 
Vb = ^V§- Ll ( 52 ) 

will remain. Therefore, to preserve the BAU, the interactions violating B and at least one of 
the Li must be out of equilibrium. In our basis-independent formulation, requiring that the 
asymmetry in B/3 — L\ be conserved will translate into a bound on the invariants 5 involving 
the Yukawa coupling A*. 

To set a bound on the L e -violating couplings or invariants, we need an estimate of the 
L e -violating rates in the early Universe at T ^ 100 GeV. If one worked in a random basis, 
one would identify /1 1 , (A^ 3 ) 1 - the 1 element of the \i and X^ q vectors, (Ag) lj : j : 1,2,3 and 
(Ag) 1J , (Ag) 1J , J : 0, ..., 3 as L e -violating couplings, and require that the rates for the processes 
mediated by these couplings to be out of equilibrium. The rate associated with any R- violating 
Yukawa coupling \r can be estimated as 

T yuk ~ 10~ 2 A^T . (53) 

and the rate for the mass corrections to gauge interaction rates can be taken as 

2 

Tmass * W- l J^g 2 T (54) 

where the inverse powers of ten account for factors of 4n, etc.. Requiring that all the "Le- 
violating" rates listed above be less than the expansion rate at T ~ 100 GeV would give the 
bounds 

\r < 10~ 7 (55) 

preserve a large enough asymmetry, since it would be proportional to (m T /T) 2 , see equations (|43| ) and (]44|). 
A similar argument can be made using the slepton mass differences in the MSSM ||, which could preserve a 
larger asymmetry if their mass differences are large and the sleptons are light enough to be still present in the 
thermal bath when the sphalerons go out of equilibrium. 

12 Note that the A^ alone, in the absence of B + L violation, is not sufficient to take the BAU to zero in the 
presence of a lepton asymmetry [|5| . This is because the lepton asymmetry carries electric charge, which must 
be cancelled by an asymmetry in the higgs and quarks. In the presence of A3, the electric charge asymmetry in 
the quarks also carries baryon number. 
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/i 1 < 100 keV (56) 

However, in the absence of a careful discussion of the basis dependence of this analysis, it is 
unclear to which elements of the coupling constant vectors and matrices ( [SUP applies, because it 
is not obvious which direction in <fi space is Li, and ([56]) is wrong because the gauge interactions 
are diagonal in any basis for <p space, and so cannot violate R parity. The correct bounds can 
be calculated by working in the T = 100 GeV mass-eigenstate basis, or by estimating the rates 
in some basis-independent way. The drawback of the thermal mass-eigenstate basis is that it is 
not the same as the zero-temperature Standard Model mass-eigenstate basis we live in today. It 
would is therefore desirable to find some basis-independent way of estimating rates associated 
with the invariants that parametize L e violation. 

We first consider which invariants must be zero for L e to be conserved. We would like to 
find a direction in <fi space that can be chosen as the left-handed lepton associated with E\. 
For this to be possible, \\ must correspond to only one plane in space, spanned by E L and 
the Higgs, and the other two A e can only intersect that plane in one direction, which is then 
the Higgs. This then means that these other two A e can also only consist of one plane each, 
because if either of these consisted of two orthogonal planes, they would intersect the A e plane 
twice. We conclude that the three A e must each consist of a single plane, and that Ag intersects 
at least one of the Ag,Ag planes, which chooses the direction for the Higgs: L\ is then the 
perpendicular direction in the Ag plane. If Ag intersects both Ag and Ag, it must do so in the 
same direction 0. The A^ and ft must have no components in the direction El, to ensure that 
L e is conserved. 

This geometry must now be translated into statements about the vanishing of certain in- 
variants. First, we need d™ ~ \bi\ 1 e \ p J 1 \ 2 = for all p, q. Recall that Ag represents a single 
plane, as in the two generation case, so is of the form (|r2|) . Setting this invariant to zero 
means that /i and X^ q do not have components along a direction in the Ag plane. Secondly, for 
similar reasons we also need 8 Pg rs = 0, V p, q, r, s, though this invariant is not independent, as 
previously discussed. Thirdly, we need 5 2 3 = for all j, to ensure that the direction in the Ag 
plane chosen to be the Higgs by \i is the same as the direction chosen by the X J e , and there 
is again a similar argument for S 3 PQ . Finally, we need 5\ 3 = 0. This implies that there is no 
Li violation among the A e : as previously argued, each A e must represent a single plane, and if 
both Ag and Ag intersect the Ag plane, it must be in the same direction. Writing Ag in the form 
(O), where H is the direction in which the three A e intersect, it is clear that AgAgAg = 0. The 
8 A J must also be zero, because the last two terms cancel the first two in the expression (|8|), in 
the absence of L\ violation. 

We now know which invariants must be zero to conserve L e exactly. However, this is not 
necessary to preserve the baryon asymmetry: we only need L e violation to be out of thermal 
equilibrium at T ~ 100 GeV. We therefore need to be able to estimate a rate for each invariant. 
It is possible to visualise any given invariant (e.g., 5i Pq ), made up of three coupling constants 
(/i, Ag, X pq in this case) as representing the squared amplitude for the process mediated by the 
"middle" interaction (Ag in this case), with thermal masses due to the "outside interactions" 
(/i, Arf in this case) on the external (j) legs i n one of the amplitudes: see figure 2. For the 



13 Note that R violation in the A e is still possible: e.g., Ag could choose the plane (0-1), Ag the plane (0-2) 
and Ag the plane (2-3). In this case, the invariants Si 23 , 6l 32 and 8 % l with i, j, k = 2,3 could still be non-zero. 
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invariant to represent a process taking place in the early Universe, the interactions appearing 
in it would have to be in thermal equilibrium. In the case of S x pq , this would require |/i| to 
be in equilibrium, so that it gives a mass to a direction in <ft space that could be the Higgs. 
Also, the row of the matrix A 1 that transforms this "Higgs" into a lepton would have to be 
in equilibrium, and we would need the component of X p ^ coupling to this lepton direction to 
be in equilibrium. We can estimate self-consistently the coupling constant for the Ag-mediated 
process to be |Ag|, as will be explained below, and the relevant component squared of X^ 9 has 
magnitude 8 pq x lA^ 7 ] 2 . We could just as well have made this argument starting with X pq , 
in which case the relevant R- violating coupling squared would be the mass 5i Pq x |/i| 2 , so the 
basis-independent bound on L e violation associated with the invariant 8i Pq is 

5\ pq x min{r, A «|,r H } < H (57) 

The rate associated with |A*| must also be in equilibrium, because otherwise no direction in </> 
space would be defined as El- 

Consider now the invariants S 2 J ' ■ These can also be imagined as corresponding to Feynman 
diagrams with thermal masses that should belong to the Higgs on both lines, and a Yukawa 
interaction that converts a Higgs into El internally. Since one of the lines is supposed to be 
El, one of the thermal masses that should belong to a Higgs but is on the El line must be out 
of equilibrium. So we obtain the similar bound 

4 i xmin{r |AJe| ,r H }</J (58) 

provided that |Ag| is in thermal equilibrium 0. 

Our previous assertion that one can estimate the coupling constant for Higgs-lepton con- 
version due to A e as |Ag| merits some further discussion. To see why this statement might 
be suspect, suppose that the Higgs is O an d E L is (pi, so that the coupling constant for the 
HElE c interaction is (Ag) 01 . If (A 1 ) 23 was much larger, it would make a greater contribution 
to |Ag|, and one could overestimate the rate for the HElE c process. However, this does not 
happen, because any element of the matrix A 1 that does not mediate the HElE c process is an 
Li-violating coupling, and is therefore required to be smaller than the Li-conserving one. 

We conclude that in basis-independent notation, the appropriate cosmological bounds are 

6l Pq xmm{T M ,T lxTl }<H (59) 

Vh r wi 



4 J xmin{r H ,r |AJ| }<if (60) 



Sl m xmm{T lx3l ,T lx7l }<H (61) 



I Ail' 



5i lK xmm{T lxil ,T im }<H (62) 

6l slpq xmm{T lX7l ,T lx7l }<H (63) 



14 It may be noted that the A e give thermal masses to both the Higgs and the leptons, so that the 4> leg with 
the thermal mass induced by A e does not have to be a Higgs, but could just as well be a lepton. This is true, 
but in the invariant 5\ l this contribution is cancelled by the second term in the invariant. So 5\ l x T\ x i\ < H 
means that the i?-violating contribution is out of equilibrium. 

16 



► X « X 



ypq ypq 



X 



Figure 2: Representation of the invariant 5^ q oc |/i^AgA^ 9 *| 2 as the Feynman supergraph for 
a scattering process mediated by A*, with thermal masses on the external legs induced by the 
interactions /x and X^ q . 

where the rates T and H are to be evaluated at T ~ T c ~ 100 GeV, using the rate estimates 
(|55|) and (0). This provides the following bounds on the invariants: 

5^|Af | 2 < 10~ 14 , p,q :1..3 (64) 

^|A^ e | 2 < 10- 14 , j:1..3 (65) 

5\ m x min{|A^ e | 2 , |A^| 2 } < 10~ 14 , j,p,q: 1-3 (66) 

8f j x min{|Ai| 2 , |A e fc | 2 } < 10~ 14 , j, k : 1..3 (67) 

C lw x min{Kf, |A^| 2 } < ID" 14 , r,s,p,q : 1..3 (68) 

where i, j, k are lepton family indices, and p, q, r, s are quark family indices. We can conserva- 
tively estimate the magnitudes of the various coupling constant vectors and matrices as taking 
their Standard Model values. If the couplings satisfy these bounds, any asymmetry in B/3 — Li 
present in thermal equilibrium in the plasma before the EPT will be preserved. Of course, it 
could just as well be an asymmetry in some other lepton flavour that is conserved, in which 
case one would obtain identical bounds, but with the lepton-flavour index "1" replaced by "2" 
or "3". 

The results ( |5"4D — (^) can be translated into basis-dependent bounds on various individual 
coupling constants, that may be more intuitive. For instance, to preserve an asymmetry in 
B/3 — Li, the R- violating coupling /i 1 should satisfy 

^h < 10- 7 (69) 

W 

in the thermal mass-eigenstate basis, where hb is the bottom Yukawa coupling. This can be 
seen by writing out the invariant 5\ 33 in the thermal mass-eigenstate basis, and requiring that 
each term in the sum making up 5\ 33 satisfy the bound (|64]). Similarly, the other L e - violating 
couplings should satisfy 

(AS 9 ) 1 < 10- 7 (70) 

(A^, (A*)' 1 < 10- 7 j:1..3 (71) 

(A 2 ' 3 ) 1J , (A 2 ' 3 ) J1 < 10- 7 J : 0..3 (72) 

in the thermal mass-eigenstate basis. 
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7 Conclusions 

In the Standard Model, there are no renormalisable lepton-number-violating operators, and 
lepton flavour number is conserved because the neutrinos are massless. However, with the 
MSSM particle content, numerous renormalisable lepton-number-violating operators can be 
constructed. This is because the leptons have the same gauge quantum numbers as a Higgs 
field, so the sleptons can interact with the same particles as the Higgs. If any of these lepton- 
number non- conserving operators are present in the Lagrangian, there is no longer a global 
symmetry defining lepton number, so a "lepton" and a "lepton-number-violating interaction" 
are no longer uniquely defined. If the baryon asymmetry we see today was present in the 
thermal soup above the electroweak phase transition, in the presence of the anomalous non- 
perturbative B + L-violating electroweak interactions, then at least one lepton flavour must 
have been conserved at the time, as we have briefly reviewed. Whether or not an asymmetry 
remains in the plasma is a physical question, which cannot depend on how one chooses to define 
a "lepton" or "lepton-number violation". 

We have therefore identified basis-independent combinations of coupling constants that 
parametrize lepton-flavour violation in the superpotential. The various couplings of the Higgs 
and the leptons can be visualised as vectors and planes in a space spanned by the Higgs and 
lepton superfields. The vectors /i and X p ^ are directions that should be the Higgs, and the \ l e 
correspond to one or two planes that should be spanned by the Higgs and the ith lepton. It 
therefore takes three couplings to construct a geometric quantity that measures R violation. 
One needs a Yukawa coupling, to provide a plane that would like to be spanned by a Higgs and 
a lepton, and then one needs two further interactions: one to choose a direction in the plane 
to be the Higgs - the perpendicular direction in the plane is now the lepton, and another to 
conflict with this choice. Our invariants correspond to the projection of a vector that should 
be the Higgs onto a direction chosen to be a lepton by other interactions. This is a geometric, 
basis-independent measure of R violation. 

We have listed all the invariants that we can construct out of three coupling constants, which 
is the minimum required, as discussed above. It is geometrically clear that these invariants are 
not all independent: indeed, there are 321 of them with three generations of leptons. Therefore, 
we have worked out how many independent invariants there are, finding the 36 expected in the 
three-generation case. In the two-generation case, we have also demonstrated explicitly how to 
express the remaining invariants in terms of the independent ones. 

In the early Universe, we can associate a lepton-flavour-violating rate with each of the 
invariants. This allows us to discuss the survival of the baryon asymmetry of the Universe in a 
basis-independent way. We have therefore determined the bounds on our invariants that would 
allow a baryon asymmetry to remain in the plasma before the electroweak phase transition. 
These bounds reduce to those of |J if flavour dependence is neglected. In the basis-dependent 
formulation that we would prefer to avoid, but which is more intuitive, we find that a baryon 
asymmetry would remain in the plasma if all the dimensionless Yukawa couplings violating 
some lepton flavour i are less than 10 -7 , and the superpotential mass term jjl 1 mixing the Higgs 
and the zth lepton satisfies hb x // < 10~ 7 /i°, where hb is the bottom yukawa, and n° mixes the 
two Higges. These are the same bounds on the Yukawa couplings as were found in previous 
work ||, K], but the bound on // is weaker. 



As noted in the introduction, we have not included the soft SUSY-breaking masses or the 
Higgs vacuum expectation value in the construction of our invariants. These would both be 
required for a complete discussion of phenomenological bounds on the i?-violating couplings in 
our basis-independent approach. It is at present unclear to us whether there is any possible 
confusion that could arise from the basis dependence of such experimental bounds, which are 
calculated in the zero-temperature mass-eigenstate basis. We expect that our invariants should 
be useful for comparing bounds on i?-violating couplings calculated at different energy scales, 
e.g., the GUT scale, the electroweak phase transition temperature, and at low energy, but such 
an analysis is left to a future paper. 

Acknowledgements 

One of us (S.D.) would like to thank Howie Haber for useful discussions. 

References 

[1] For reviews, see, for instance: 

H.P. Nilles, Phys. Rep. 110 (1984) 1 ; 

H.E. Haber and G.L. Kane, Phys. Rep. 117 (1985) 75 ; 

G.G. Ross, Grand Unified Theories, (Benjamin-Cummings, Menlo Park, CA, 1985); 

R. Barbieri, Riv. Nuovo Cimento 11, 1 (1988). 

[2] G. Farrar and P Fayet, Phys. Lett. B76 (1978) 575; 
S. Weinberg, Phys. Rev. D26 (1982) 287; 
N. Sakai and T. Yanagida, Nucl. Phys. B197 (1982) 133 . 

[3] See, for instance, the following papers and references therein: 

G . Bhattacharyya, D. Choudhury and K. Sridhar, Mod. Phys. Lett. A10 (1995) 1583; 

G . Bhattacharyya, D. Choudhury and K. Sridhar, Phys. Lett. B355 (1995) 193; 

C.E. Carlson, P. Roy and M. Sher, Phys. Lett. B357 (1995) 99; 

M. Hirsch, H.V. Klapdor-Kleingrothaus and S.G. Kovalenko, Phys. Rev. Lett. 75 

(1995) 17 ; 

K.S. Babu and R.N. Mohapatra, Phys. Rev. Lett. 75 (1995) 2276; 

A. Yu. Smirnov and F. Vissani, IC-96-16, |hep-ph/9601387| ; 



D. Choudhury and P. Roy, Phys. Lett. B378 (1996) 153 

M. Chaichian and K. Huitu, HU-SEFT-R-1996-09 and |hep-ph/96034T2 



H. Nunokawa, A. Rossi and J.W.F. Valle, FTUV-96-33, |hep-ph/9606445| . 

[4] BA. Campbell, S.Davidson, J. Ellis and K.A. Olive, Phys. Lett. B 256 (1991) 457, 
Astroparticle Phys. 1 (1992) 77. 

[5] H. Dreiner and G.G. Ross, Nucl. Phys. B410 (1993) 188. 



[6] S. Davidson and J. Ellis, CERN-TH/96-258, CfPA-96-th-20, MPI-PhT-96-68, hep- 
ph/9609451 



19 



[7] G. 't Hooft, Phys. Rev. Lett. 37 (1976) 8; 

F.R. Klinkhamer and N.S. Manton, Phys. Rev. D30 (1984) 2212; 
P. Arnold and L. McLerran, Phys. Rev. D36 (1987) 581. 

[8] J.R. Espinosa, Nucl. Phys. B475 (1996) 273; 

M. Carena, M. Quiros and C.E.M. Wagner, CERN-TH-96-30, |hep-ph/9603420| ; 



D. Delepine, J.M. Gerard, R. Gonzalez Felipe and J. Weyers, UCL-IPT-96-05, |hep- 
ph/ 96044401 ; 



J. Cline and K. Kainulainen, |hep-ph/9605235 . 



[9] S.Dimopoulos and L.J. Hall, Phys. Lett. B 196 (1987) 135; 
J. Cline and S. Raby, Phys. Rev. D 43 (1991) 1781; 

R. J. Scherrer, J. Cline, S. Raby and D. Seckel, Phys. Rev. D44 (1991) 3760; 
A. Masiero and A. Riotto, Phys. Lett. B289 (1992) 73. 

[10] V.A. Kuzmin, V.A. Rubakov and M. E. Shaposhnikov, Phys. Lett. 155B (1985) 36. 

[11] S. Y. Khlebnikov and M.E. Shaposhnikov, Nucl. Phys. B308 (1988) 885; and Phys. 
Lett. B387 (1996) 817. 



[12] S. Davidson and R. Hempfhng, MPI-PhT-96-89, |hep-ph 9609497 



[13] B.A. Campbell, S.Davidson, J. Ellis and K.A. Olive, Phys. Lett. B 297 (1992) 118; 
J. M. Cline, K. Kainulainen and K. A. Olive, Phys. Rev. Lett. 71 (1993) 2372, Phys.Rev. 
D49 (1994) 6394. 

[14] J. Harvey and M. Turner, Phys. Rev. D42 (1990) 3344. 

[15] L.D. Landau and E.M. Lifshitz, Statistical Physics (Pergamon Press, 1989). 

[16] S. Davidson, K. Kainulainen and K.A. Olive, Phys. Lett. B335 (1994) 339. 



20 



